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Abstract 

We report the first measurement at the LHC of coherent photoproduction of p° mesons in ultra¬ 
peripheral Pb-Pb collisions. The invariant mass and transverse momentum distributions for p° pro¬ 
duction are studied in the 7T + 7T~ decay channel at mid-rapidity. The production cross section in the 
rapidity range |y| < 0.5 is found to be da/Ay = 425 ± 10 (stat.) (sys.) mb. Coherent p° produc¬ 

tion is studied with and without requirement of nuclear breakup, and the fractional yields for various 
breakup scenarios are presented. The results are compared with those from lower energies and with 
model predictions. 
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1 Introduction 

Charged particle beams at the LHC generate an electromagnetic field which can be regarded as a beam 
of quasi-real photons; thus at the LHC, besides hadronic interactions, also photonuclear and photon- 
photon interactions occur. Collisions in which the impact parameter exceeds the sum of the radii of 
the incoming beam particles are called ultra-peripheral collisions (UPC). In UPC the cross section for 
hadronic processes is strongly suppressed, while the cross sections for two-photon and photonuclear 
interactions remain large. This is particularly the case for heavy ions, because the intensity of the photon 
flux grows with the square of the ion charge, Z. A number of reviews of UPC are available; e.g., [1, 2]. 
The ALICE Collaboration has previously studied exclusive photoproduction of 7/t/r in ultra-peripheral 
Pb-Pb and p-Pb collisions [3-5]. 

Exclusive photoproduction of p° vector mesons, Pb + Pb —> Pb + Pb + p°, can be described as the fluc¬ 
tuation of a quasi-real photon into a virtual p° vector meson, which then scatters elastically off the target 
nucleus. Two cases can be distinguished. When the interaction involves the complete target nucleus, the 
process is called coherent. In this case the target nucleus normally remains intact. If the virtual p° vector 
meson scatters off only one of the nucleons in the target, then the process is called incoherent and in this 
case the target nucleus normally breaks up, emitting neutrons at very forward rapidities. For coherent 
processes, the size of the lead ion restricts the mean transverse momentum of the vector meson to be 
about 60 MeV/c corresponding to a de Broglie wavelength of the nuclear size, while it is of the order of 
500 MeV/c for incoherent processes. 

Because of the strong electromagnetic fields in ultra-peripheral collisions of heavy ions, multiple pho¬ 
tons may be exchanged in a single event. The additional photons can lead to excitation of the nuclei. 
The dominant process is the excitation to a Giant Dipole Resonance [ 6 ]. As these photonuclear pro¬ 
cesses occur on a different time scale, they are assumed to be independent, so the probabilities factorize. 
The excited nucleus typically decays by the emission of neutrons at very forward rapidities. The sig¬ 
nature of these processes is thus a p° vector meson with very low transverse momentum which may be 
accompanied by a few neutrons at very forward rapidities but no other particles. 

Photoproduction of p° vector mesons on nuclear targets has been studied in fixed target experiments with 
lepton beams [7], and more recently in ultra-peripheral collisions by the STAR Collaboration at RHIC 
at yAnn = 62 [ 8 ], 130 [9], and 200 GeV [10]. STAR has also observed coherent photoproduction of the 
p°(1700) [11]. 

The p° vector meson gives the dominant contribution to the hadronic structure of the photon. For proton 
targets, the process y+ p —»• p° + p contributes about 10 % to 20 % of the total y+p cross section, 
depending on energy [12]. Scaling from a nucleon target to a nuclear target is often done using the 
Glauber model assuming Vector Meson Dominance [13]. The large value of <7(7 + p — > p° + p) means 
that for heavy nuclei one may reach the limit where the target appears like a black disk and the total 
p° + A cross section approaches 2kR 2 a ( Ra is the nuclear radius). The situation may, however, be more 
complicated for several reasons. The cross section <7(7 +p — > p° + p) has contributions both from 
Reggeon and Pomeron exchange, and its energy dependence is therefore not monotonic. Furthermore, 
the nuclear medium might modify the Reggeon and Pomeron components differently. There may also be 
interference between the p and p' production amplitudes, and these amplitudes may be affected by the 
nuclear environment in a different way [14]. A detailed discussion of models for photoproduction of p° 
on complex nuclei based on data from fixed target experiments can be found in [13]. 

The cross sections measured by STAR [8-10] at RHIC were found to be about a factor two less than that 
predicted by the calculation of Ref. [15], while in agreement with STARLIGHT [16]. The reason for the 
difference between these two models, which both use the Glauber model to obtain the 7 -nucleus cross 
section, will be discussed below. The many issues associated with calculating the photonuclear p° cross 
section and the discrepancies between models thus call for more data. In particular, it is important to 
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establish if the trends seen at lower energies persist at higher energies. 

Moreover, the total cross section for exclusive p° production is very large at LHC energies, with the 
models mentioned above predicting that it could be between 50-100% of the total hadronic inelastic 
cross section. It could thus constitute a significant background, e.g. at the trigger level, to low multiplicity 
peripheral hadronic interactions and to other types of ultra-peripheral collisions. It therefore has to be 
well understood. The high statistics in the p° sample allows the predictions for exclusive p° production 
accompanied by nuclear fragmentation to be tested with good precision. 

This paper presents the first measurement of the cross section for coherent photoproduction of p° vector 
mesons in Pb-Pb collisions at the LHC. The p° is reconstructed using the n n decay channel in the 
rapidity range |y| < 0.5. The rapidity interval corresponds to a y-nucleon center of mass energy in 
the range 36 < VV/v < 59 GeV with (VTy ; y) = 48 GeV, about a factor of 4 higher than in any previous 
measurement [10]. The cross section is measured for the cases of no neutron emission and for at least 
one emitted neutron. The new data presented in this paper will hopefully help to clarify some of the 
theoretical uncertainties mentioned above. 

2 The ALICE experiment and the UPC trigger 

A full description of the ALICE detectors and their performance can be found in [17, 18]; here, only the 
components relevant for this analysis will be briefly described. The Inner Tracking System (ITS) and 
Time Projection Chamber (TPC) are used to measure and identify the tracks of the decay products of the 
p° vector meson. The ITS consists of six layers of silicon detectors covering the full azimuthal angle. 
The two innermost layers form the Silicon Pixel Detector (SPD) with a pseudorapidity acceptance of 
1 77 1 < 1.4. The SPD also provides trigger information at the lowest level. Two layers of silicon drift and 
two of silicon strip detectors complement the ITS, and all six layers have an acceptance of 1 77 1 < 0.9. 
The TPC is the main tracking detector of ALICE. It has a Ne-C02-N2 gas mixture contained in a large - 
almost 90 m 3 - cylindrical drift detector with a central membrane at high voltage and two readout planes, 
composed of multi-wire proportional chambers, at the end caps. It covers the full azimuth and 1 77 1 < 0.9 
for full length tracks. It also provides a measurement of the ionization energy loss, d E /dr, which allows 
the identification of particles. The TPC and ITS are situated inside a large solenoid magnet providing a 
B = 0.5 T field. 

The measurement of neutrons emitted at forward rapidities is performed with a set of two neutron Zero 
Degree Calorimeters (ZDC) located 114 m away on each side of the interaction point. The ZDC has 
a 99% detection probability for neutrons with 1 77 1 > 8.8[19]. Figure 1 illustrates the capabilities of the 
ZDC to separate the emission of zero, one or several neutrons at zero degrees. The sample appearing in 
this figure was obtained from events fulfilling the event selection described in Section 3. 

In addition to the SPD, this analysis uses the Time of Flight (TOF) and VZERO detectors for triggering. 
TOF is a large cylindrical barrel surrounding the TPC. It has 18 sectors in azimuth, each made of multigap 
resistive plate chambers distributed in five gas-tight modules, with a total of 152928 read-out channels 
and an intrinsic time resolution better than 50 ps. The pseudorapidity acceptance is the same as for the 
TPC. The VZERO consists of two arrays of scintillators called VZERO-A and VZERO-C, covering the 
pseudorapidity ranges 2.8 < r\ < 5.1 (VZERO-A) and —3.7 < 17 < —1.7 (VZERO-C). Its time resolution, 
better than 500 ps [20], allows beam-beam collisions to be distinguished from beam-gas collisions. 

The data used for this analysis were collected during the 2010 Pb-Pb run of the LHC at an energy of 
y(?NN = 2.76 TeV. Two different triggers were used. At the beginning of the run, when the luminosity 
was low, the trigger requirement was at least two hits in the TOF detector. When the luminosity was 
increased the trigger selection was strengthened to improve the purity by additionally requiring at least 
two hits in the outer layer of the SPD, and no activity in any of the VZERO arrays. 
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Fig. 1: Energy deposit in the Zero-Degree Calorimeters. The curves correspond to Gaussian fits for 0, 1,2, 3, or 4 
neutrons entering the calorimeter. The plot is for events satisfying the requirements (i)-(vi) described in the text. 

The luminosity is determined from the cross section for triggering on at least one neutron in the ZDC de¬ 
tectors [19]. This cross section has been determined from a van der Meer scan [21] to be 
371.4 ± 0.6 (stat.)^ (syst.) b [18]. The integrated luminosities for the two samples are 48 1 2 mb -1 
(TOF trigger only) and 214/J/ mb 1 (SPD+TOF+VZERO trigger). 

3 Track and event selection 

In addition to the trigger selection, the events used for the analysis are required to fulfill the following 
requirements: 

i) a primary vertex has to be identified within 10 cm of the nominal interaction point position, along 
the beam direction; 

ii) the event is required to have exactly two tracks reconstructed in the TPC and ITS satisfying the track 
selections discussed below; 

iii) the VZERO arrays are required to be empty (the difference between the offline and online VZERO 
selection will be discussed below); 

iv) the energy loss in the TPC has to be consistent with that for pions within 4 standard deviations from 
the Bethe-Bloch expectations, i.e., Ao 2 + A a 2 - < 16 (see Fig. 2); 

v) the track pairs used to define the coherent signal have to have a transverse momentum below 
150 MeV/c and rapidity |y| < 0.5, the latter requirement being imposed to avoid edge effects; 

vi) the track pairs used to define the coherent signal are required to have tracks of opposite charge. 

The background estimated from like-sign pairs {n + n + and 71 n ) is below 2% and it is subtracted from 
the final sample bin-by-bin in invariant mass. 

The track selection requires that each track has at least 70 space points, out of a maximum of 159, 
in the TPC and a % 2 per degree of freedom from the Kalman fit procedure better than 4. Each track 
has at least one hit in the SPD with a % 2 per ITS hit less than 36. The distance of closest approach 
between the track and the primary vertex has to be less than 2 cm along the beam direction and less than 
0.0182 +0.035/p}- 01 cm (pj in GeV/c) in the plane perpendicular to the beam direction. These track 
selection cuts are based on studies of the detector performance [18]. 
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Three other track selections are used in order to estimate systematic errors. These differ from the default 
track selection described above in the following ways: (a) accepting tracks reconstructed only in the 
ITS in addition to combined ITS-TPC tracks satisfying the default track selection; (b) using only TPC 
information and accepting tracks having at least 50 space points in the TPC; (c) using the default track 
cuts with stronger requirements on TPC variables. The latter requirements meant that the tracks had to 
pass at least 120 of the 159 TPC pad rows and have a cluster in more than 80% of the crossed pad rows. 
For the cross section calculation, the mean of the results of the four different track selection methods is 
used. The systematic error related to the track selection is estimated from the deviation from the mean. 
This contributes + 3 7 0 % to the systematic error. 

The momentum resolution of the ALICE central barrel tracking system [18] translates into a resolution 
in transverse momentum of single Ji 1 71 -pairs better than 4 MeV/c in the kinematic range studied here. 
Similarly, the resolution in invariant mass varies between 2 MeV/c 2 (M nn = 0.4 GeV/c 2 ) and 6 MeV/c 2 
(M xk = 1.5 GeV/c 2 ). 
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Fig. 2: Identification of K 1 from the TPC d/i/cL. The x- and y-axes show the deviation of the measured energy 
loss from the Bethe-Bloch (BB) expectations for positive and negative tracks, respectively. The scale is such that 
one unit corresponds to one standard deviation. The circle, corresponding to 4(7, shows the selection used. The 
entries in the upper right corner are from e + e~ pairs produced in two-photon interactions. The plot is for events 
satisfying the requirements (i)-(iii) and (v) - (vi) described in the text. 

The ionization energy loss for the selected tracks is shown in Fig. 2. The scale on both axes is in units 
of the number of standard deviations from the Bethe-Bloch expectation in the TPC; in this way the 
dependence on track momentum is removed. Pions can be clearly identified by the 4 <7 circle centered 
on ( 0 , 0 ), while the events above and to the right of the pions are mostly e e pairs from yy —> e + e . 
This figure shows that any possible contamination from kaons or protons in the sample is negligible. 
There could, however, be a contamination from muons from the process yy—^ /.( 1 /J , which cannot be 
distinguished from pions using the energy loss. This contribution can be estimated from the number of 
e + e~ pairs in the data sample, as the cross sections for yy —^ /r /i and yy —>• e ] e are about the same at 
midrapidity for invariant masses well above threshold. It can also be calculated from STARLIGHT [22, 
23]. Both methods give an expected number of muon pairs of about 5%, which is not corrected for, but 
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added to the systematic error. The contribution from yy—>• Jl 1 71 is expected to be much smaller than 
from yy —> /i 1 /i . The K + n~ cross section is reduced by the form factor of the pion, see e.g. [24], so 
this contribution is not considered. 

4 Data analysis 

Using the event and track selection described in the previous section, the four-momenta of the two tracks 
are constructed and pair variables are extracted. The resulting distribution of the pair transverse momen¬ 
tum is shown in Fig. 3 for events with 0.4 < M nn <1.1 GeV/c 2 and |y| < 0.5. A peak at low transverse 
momentum (pj < 0.15 GeV/c), corresponding to coherent production, is clearly seen. The distribution 
is compared with the corresponding distributions from STARLIGHT [16, 23] events for coherent and 
incoherent p° production, processed through the detector response simulation based on GEANT 3. The 
coherent peak is shifted to slightly lower pj in data than that predicted by STARLIGHT. A similar trend 
has been observed by STAR at lower energies [25]. The shape of the coherent peak in the p-\ distribution 
is determined by the nuclear form factor. The form factor used in STARLIGHT is consistent with what is 
obtained from elastic electron-nucleus scattering, which probes the charge content of the nucleus. Since 
the p° couples to both neutrons and protons, a possible explanation of this difference could thus be the 
presence of a “neutron skin”. The effect, however, appears larger than what the current limit on the dif¬ 
ference between neutron and proton radius in 208 Pb (0.3 fm) allows [26], and is thus not fully understood. 
Data also show a dip around pj = 0.12 GeV/c, which is not present in the model. The absence of this 
dip in the model can be understood from the fact that in STARLIGHT the transverse momentum of the 
photon is considered, and this reduces the dip one would expect from the form factor of the target nu¬ 
cleus alone. In a Glauber calculation, the transverse momentum distribution is determined from a Fourier 
transform of the nuclear profile function, see e.g. [2], and the direct dependence on the form factor is 
only an approximation; this could also contribute to explaining the difference between STARLIGHT and 
data. The high-pj tail of the distribution is very well described by the incoherent pj spectrum from 
STARLIGHT. 

The transverse momentum distribution for coherent production may also be parameterized as an expo¬ 
nential, dN/dt o< exp (At) where t = —pj. Fitting the ALICE data to such a function gives b = 428 ± 
6(stat.) ± 15(syst.) GeV 2 /c 2 . The systematic error has been obtained as the difference in slope between 
STARLIGHT events and STARLIGHT events processed through the full detector simulation. The AL¬ 
ICE result can be compared with the corresponding measurement by STAR, where b = 388 ± 24 GeV 2 /r 2 
was found [10]. The STAR and ALICE results are consistent within errors if one takes into consideration 
that b is expected to be ~4-8% larger for a lead nucleus than for a gold nucleus because of the difference 
in size (one expects b <=< R 2 ). The fit was performed for \t\ > 0.002 GeV 2 /c 2 to avoid interference effects 
at very low pr [10]. 

The final sample of coherent p° —>■ n + n~ candidates is corrected for acceptance and efficiency in invari¬ 
ant mass bins. The event sample used to determine the corrections has uniform distributions in invariant 
mass, rapidity, transverse momentum, and azimuthal angle over the ranges 2 m n < M nn <1.5 GeV/c 2 , 
|y| < 1.0, pt < 0.15 GeV/c, and 0 < (j) < 2u. Using a flat distribution in transverse momentum is 
justified over the narrow range p-\ < 0.15 GeV/c, where the acceptance and efficiency are constant. All 
models predict only a very small variation of the cross section over the range |y| <0.5 (see Fig. 5 below) 
so also for rapidity a uniform input distribution is justified. The advantage of using a flat input distribu¬ 
tion in invariant mass is to obtain sufficient statistics in the tails of the distribution. If one were to use a 
p°-shape as input, one would need enormous statistics to cover the high and low invariant mass ranges. 
The p° candidates are assumed to be transversely polarized. This is expected from helicity conservation 
and has been confirmed by photoproduction measurements [10, 27]. This polarization translates into a 
dn^/dO, °< sin 2 (0) angular distribution of the K ' n decay products in their center of mass system (0 is 
here measured relative to the direction of flight of the p° in the y-nucleon center of mass system). All 
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Fig. 3: Transverse momentum distributions for K 1 n -pairs. The dashed (red) and dash-dotted (blue) histograms 
show the normalized pj distribution from STARLIGHT passed through the detector response simulation for co¬ 
herent and incoherent p° production, respectively. The solid (black) histogram is the sum of the two. 

generated samples serve as input to a full detector simulation using GEANT 3 for the propagation of 
particles through the detector. Selection criteria are applied in the same way as done for real events. The 
variation of the detector configuration during the data taking period is included in the detector response 
simulations. The product of acceptance and efficiency varies from about 2% at the low end of the studied 
invariant mass interval ( M nn = 0.6 GeV/c 2 ) to about 12% at the high end (M nn = 1.5 GeV/c 2 ). 



Fig. 4: Invariant mass distribution for pions corrected for acceptance and efficiency. The solid (blue) curve corre¬ 
sponds to a fit to the Soding parameterization Eq. 1, and the dashed (green) curve shows the resonant contribution 
only. The dot-dashed (red) curve shows the fit to the Ross-Stodolsky parameterization Eq. 3. The parameters of 
the fit are given in the text. 

The uncertainty in the trigger efficiency is obtained by comparing the measured trigger efficiency with 
the one in the detector response simulation in a data sample taken with a ZDC trigger [4]. The result 
is a trigger efficiency uncertainty of In addition, a correction is applied for the trigger dead time 

resulting from after pulses in the TOF. The systematic error on this correction is estimated to be ±1.3%. 

The invariant mass distribution of the p° candidates, corrected for acceptance and efficiency and nor¬ 
malized by the luminosity to provide a cross section, is shown in Fig. 4. It is well known that the shape 
of the p° in photoproduction deviates from a pure Breit-Wigner resonance [7-10, 27]. Several different 
parameterizations exist to describe the shape, with one of the most often used being a formula due to 
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Soding, where a continuum amplitude, B , is added to the Breit-Wigner resonance [28]: 


da 


M 7C7t M p oT(M' 


nn J 


Ml n -Ml 0 + iM p0 T(M nn ) 


+ B 


( 1 ) 


Here, A is the amplitude of the Breit-Wigner function, B is the amplitude of the direct non-resonant 
7r + 7T~ production, and the mass dependent width is given by 


M n o 

r(M^ w ) = r p o— 

ivljin 

with m K the mass of the pion. Eq. 1 was fitted to the measured M nn distribution with M p o, I ' p a, 
A, and B as free parameters. The fit gives M p o = 761.6 ± 2.3 (stat.)^^ (syst.) MeV/c 2 and I ' p n = 
150.2 ± 5.5 (stat.)])^ 0 (syst.) MeV/c 2 , in agreement with the values reported by the PDG [29]. The 
ratio of the non-resonant and resonant amplitudes is found to be \B/A\ = 0.50 ± 0.04 (stat.)^'^ (syst.) 
(GeV/c 2 ) -1 / 2 . The systematic errors are obtained by varying the fitting method (j} or log likelihood 
minimization), track selection (as discussed above), and fitting range. 

The ratio \B/A\ is lower than what was found by STAR with Au targets, \B/A\ = 0.81-0.89 (GcV/c 2 ) '/ 2 
for (Wyff) in the range 7-12 GeV [8-10]. The result from ZEUS with proton targets for (Wy^) in the 
range 55-90 GeV shows that \B/A\ varies with the y-proton momentum tranfer [27]. The average is 
\B/A\ = 0.67 ±0.02 (stat.)±0.04 (syst.), while for momentum transfers of the same order as for coherent 
production \B/A\ ~ 0.8. The lower value of \B/A\ observed by ALICE may indicate that the non¬ 
resonant contribution is more strongly absorbed in heavy nuclei at high energies, as had been previously 
suggested [30]. 

Other parameterizations of the p° shape are possible, and as a cross check the invariant mass distribution 
was also fit to a Ross-Stodolsky function [27, 31]: 


44 nn 4m n 

M 2 - 4m 2 


da 

dM KK 


f 


\J M K7l M p oT(M 

11 


— M 2 p0 ± iM p aT{M n 



with a slightly different definition of the mass dependent width 


r(M**) — TpO 


M 1 

lrt TT 


M 2 0 
L P° 


4m l 

4 ml 


( 3 ) 


( 4 ) 


As can be seen in Fig. 4, this parameterization also described the observed shape of the invariant mass 
distribution well and gave a p° mass (M p o = 769.2 ± 2.8 (stat.)^ ® (syst.) MeV/c 2 ) and width (Fpo = 
156.9d= 6.1 (stat.)^ 5 7 9 3 (syst.) MeV/c 2 ) consistent with the PDG values. The deviation from apure Breit- 
Wigner shape is given by the parameter k, which was found to be k = 4.7 ±0.2 (stat.)] [[ 9 (syst.). This can 
be compared to the corresponding value for proton targets from ZEUS [27] and HI [32] at HERA. ZEUS 
finds k = 5.13 ±0.13 averaged over all momentum transfers and k ~ 6 for t = 0, while HI reports k = 
6.84 ± 1.00 averaged over all momentum transfers. The larger value of k for proton targets again indicates 
that the invariant mass distribution for Pb-targets deviates less from a pure Breit-Wigner resonance, as 
was also found using the Soding formula. 

As can be seen in the lower part of Fig. 4, there is a hint of a resonance around 1.3 GeV/c 2 . This 
may be understood from two-photon production of the /2(1270) meson followed by its decay into two 
pions, y±y->/2(1270) —>• K + 7i~. This meson is a “standard candle” in two-photon interactions with a 
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well known 77 coupling, but it has so far not been observed in ultra-peripheral collisions because of the 
large background from photonuclear processes. The significance of the excess over the p° Breit-Wigner 
distribution is estimated to be 4 1 2 , where the error comes from the uncertainty in the skewness of the 
Breit-Wigner distribution (parameter k in the Ross-Stodolsky formula). 

The normalized yield of p°s (/V y j e id) is obtained by integrating the resonant part of Eq. 1 (obtained by 
setting B = 0 and taking the other parameters from the fit) from 2 m n to 1.5 GeV/c 2 . The systematic 
error on the number of extracted p°s is obtained by varying the fitting method (% 2 or log likelihood 
minimization) and fitting range, resulting in an error of /j 4 %. The uncertainty in the track selection 
gives an additional error of +3 q% as discussed above. Both Eq. 1 and 3 describe the observed shape 
equally well (the integrated yield differ by less than 0.5%), so no additional systematic error was added 
to the yield because of the choice of fitting function. 

It is worth noting that the shape of the resonant contribution (shown by the dashed curve in Fig. 4) is quite 
different from the shape of the measured n + n~ invariant mass distribution. However, the integrated yield 
between 2m n and 1.5 GeV/c 2 does not deviate by more than around 1% if the non-resonant amplitude is 
included in the integration. 

The number of extracted p°s is corrected for the following 3 contributions: incoherent events with 
p-\ < 0.15 GeV/c ( /incoh), events which have one or more additional SPD tracklets (/spd), and the number 
of coherent p° events lost by the VZERO offline timing requirement (/v zero)- 

The number of incoherent events with pj < 0.15 GeV/c is estimated in two different ways: first fitting 
the sum of two exponentials in p\ to the pj distributions and integrating the fitted functions over the 
interval chosen for the coherent selection (pj <0.15 GeV/c), and second using the fit to the STARLIGHT 
templates shown in Fig. 3. The correction for incoherent events is found to be 5.1% in both cases with 
an uncertainty estimated from using different track selections of ±0.7%. 

The track selection (a) above allows one to check the events for any additional activity in the ITS, for 
example from tracks with low momenta, which do not reach the TPC, using SPD tracklets, defined as 
any combination of hits from the two SPD layers. Rejecting events with one or more extra tracklets, not 
associated with the two good tracks coming from the primary vertex, removes 3.0% of the events in the 
signal region. Since true UPC events should have no additional tracks, the extracted yield is corrected 
for this. In the Monte Carlo samples of coherently produced p°s, the same cut removes only 0.5% of the 
events which is taken as the systematic error associated with this cut. 

The events selected by the SPD+TOF+VZERO trigger are required to have no online signal in the 
VZERO detector. A similar cut is also applied offline to the events triggered by TOF only. The VZERO 
offline selection is further refined using the timing information. This selection has been tuned to work 
well for hadronic interactions, which typically have a non-zero signal in the VZERO on both sides. In 
the ultra-peripheral events studied here, where the VZERO is required to be empty, the offline selection 
is less reliable, and a coherent signal can be observed in the events with 2 tracks rejected by the offline 
VZERO requirement. The increase in the coherent signal when the offline VZERO selection is not used 
amounts to 10.0%. The systematic error of this number is obtained from the estimated contamination 
from hadronic events following from this looser cut. This contamination is determined from the fraction 
of the events which have a signal in the ZDCs, resulting in a systematic error of 

The corrected number of coherent p°s is then obtained from 

M coh = _ N ? ield _ (S) 

P 1 + /incoh ± /SPD ± fw ZERO ’ 

with /ncoh = 0.051 ±0.007, /spd = 0.030±0.005, and /v zero = — 0.100/qqqq. From this number the 
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Variable 

Systematic error 

Luminosity 

+6.5% 

-5.1% 

Trigger efficiency 

+3.8% 

-9.0% 

Trigger dead time correction 

±1.3% 

Signal extraction 

+0.8% 

-1.4% 

Track selection 

+3.7% 

-3.0% 

Particle ID 

+0.0% 

-5.0% 

Incoherent contribution 

±0.7% 

SPD tracklets 

±0.5% 

VZERO offline selection 

+0.0% 

-3.1% 

Total 

+9.2% 

-11.2% 


Table 1: Summary of the systematic error in the cross section calculation. The numbers are for the 
SPD+TOF+VZERO trigger sample. For a discussion of the TOF only trigger sample and the separation between 
correlated and uncorrelated errors of the two samples, see the text. 

differential cross section is calculated as 

— = -. ( 6 ) 

dv L jnt • Ay 

The systematic errors discussed above are summarized in Table 1. They have been evaluated for the 
SPD+TOF+VZERO trigger sample, which contains more than 80% of the total integrated luminosity. 
The total error is obtained by adding the individual errors following the description in [33]. The two 
trigger samples, with appropriate errors, are compared as a cross check. They use different trigger com¬ 
binations and were taken under quite different running conditions, with the typical hadronic minimum 
bias interaction rate being around 10 Hz during the early part of the run when the TOF only trigger was 
used and around 200 Hz during the later part of the run when the SPD+TOF+VZERO trigger was used. 
The correction factor for trigger dead time due to after pulses was thus very different for the two samples 
(~ I during the early part and ~5 during the later part). 

To make a comparison of the cross sections measured under the different trigger conditions, the system¬ 
atic eiTors are separated into correlated and uncorrelated errors for the two trigger samples. The fully 
correlated errors are those related to luminosity, incoherent contribution, trigger efficiency, and particle 
identification. The fully uncorrelated errors are those related to the VZERO offline selection (different 
VZERO thresholds were used for the two data samples), the cut on SPD tracklets, and trigger dead time. 
The errors related to the signal extraction and track selection are found to be partly correlated, but are 
decorrelated for the comparison. This gives a cross section da/dv = 466 ^ mb for the sample taken 
with the TOF only trigger and da/dv = 4I4 1 ]/ mb for the sample taken with the SPD+TOF+VZERO 
trigger. The error is obtained from the squared sum of the statistical and uncorrelated systematic error. 
The difference of 12% corresponds to 1.8 standard deviations. The final cross section is obtained as 
the weighted mean of the cross sections of the two samples. The weighting procedure provides a total 
eiTor, including both the statistical and uncorrelated systematic components. The uncorrelated compo¬ 
nent is separated from the total error by subtracting in quadrature the error obtained in the case when 
only the statistical errors are used for the weighting. The uncorrelated systematic error is then added 
in quadrature to the correlated systematic error to obtain the total systematic error. The final result is 
da/dy = 425 ± 10(stat.);/jQ (syst.) mb. 
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Fig. 5: The cross section for coherent photoproduction of p°, dc/dy, in ultra-peripheral collisions for the three 
models compared with the ALICE result. 

In addition to the p° cross section, the cross section for two-photon production of e + e~ pairs in the 
range 0.6 < M ee < 2.0 GeV/c 2 and 1 171 2 1 < 0.9 ( 171.2 are the pseudorapidities of the two tracks) was 
measured. The analysis is similar to the one for p° but the PID requirement was modified to accept 
electrons rather than pions. The detector efficiency is determined using STARLIGHT events processed 
through the full ALICE detector simulation. The result is a(0.6 < M ee < 2.0 GeV, 1 771 n 2 1 < 0.9) = 9.8 ± 
0.6 (stat.)+j T (syst.) mb, which is in good agreement with the STARLIGHT [22] prediction for the same 
selection in invariant mass and pseudorapidity (a = 9.7 mb). The cross sections for the individual 
trigger samples are 11.8± 1.6 (stat.)"]^ (syst.) mb (TOF only trigger) and 9.4±0.7 (stat.)^^ (syst.) mb 
(SPD+TOF+VZERO trigger). 

As discussed above, photoproduction of vector mesons may occur in interactions where additional pho¬ 
tons are exchanged between the nuclei, leading to neutron emission in the forward region. These neutrons 
may be detected in the ALICE ZDCs. Four Gaussian distributions centered around each peak with means 
and variances constrained to x n = nx\ and o n = y/nC\ have been fitted to the ZDC energy distribution 
shown in Fig. 1 . Here, x\ and Ci are the position and width of the peak corresponding to one neutron, 
and n is the number of neutrons. In order to separate different cases of neutron emission, the minima 
between the first three Gaussians are used. The minimum between zero and one-neutron emission lies at 
half the energy per nucleon and it is roughly three sigma away from the adjacent peaks. A given event is 
considered to have no neutron in the ZDC if the energy registered in the calorimeter is less than 600 GeV, 
one neutron if the energy lies between 600 GeV and 2000 GeV and more than one neutron if the energy 
is above 2000 GeV. 

The events are divided into different groups as follows: no neutrons emitted in any direction (OnOn), 
at least one neutron emitted in any direction (Xn), at least one neutron emitted in one direction and no 
neutron emitted in the other direction (OnXn), at least one neutron emitted in both directions (XnXn). 

The corrections applied in obtaining the cross section from the measured yield are independent of the 
ZDC signal. The fractional yield for each fragmentation selection thus reflects the relative p° production 
cross section. The only exception to this is the correction for the incoherent contribution (/i nco h), which 
is expected to be higher when a signal is required in the ZDCs. This correction is thus calculated for 
each ZDC selection separately, using the same method as described above. 
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Fig. 6: Excitation function for coherent and exclusive p° production. The results from ALICE and STAR[8-10] 
are compared with the STARLIGHT and GDL predictions for Pb-Pb and Au-Au. 

5 Results and discussion 

The coherent p° photoproduction cross section, da/dy, is shown and compared with model predictions 
in Fig. 5. The measured cross section is in agreement with STARLIGHT [16] and the calculation by 
Gon 9 alves and Machado (GM) [34], while the GDL (Glauber-Donnachie-Landshoff) prediction [15, 35] 
is about a factor of 2 higher than data. The calculation by GM is based on the Color Dipole model, while 
STARLIGHT and GDL use the photon-proton cross section a(y + p —> p° + p) constrained from data 
as input. In STARLIGHT, the y-nucleon cross section is given by the parameterization a = 5.0 IT^ 20 + 
26.0 W^ 1 ' 23 pb (WyN in GeV), while GDL use the Donnachie-Landshoff model [36] for the total pN cross 
section. All calculations use the Glauber model to scale the cross section from y-nucleon to y-nucleus. 

The STAR Collaboration has published the total coherent p° photoproduction cross section at three 
different energies [8-10]. To be able to compare the current result to those, one has to integrate da/dy 
over the whole phase space, which can only be done using models. The extrapolation factor from |y| < 
0.5 to all rapidities is calculated as the mean of the values obtained from the STARLIGHT (10.6) and GM 
(9.1) models, and the deviation of the two from the mean (« 8%) is added to the systematic error. This 
gives a(Pb + Pb —y Pb + Pb + p°) = 4.2±0.1(stat.)^Qg(syst.) b at ^/.vnn = 2.76 TeV. The total cross 
section as a function of ^Ann is shown in Fig. 6, where the results from AFICE and STAR Collaborations 
are compared with the STARFIGHT and GDL calculations. The total cross section increases by about a 
factor of 5 between the top RHIC energy and ^sxn = 2.76 TeV. 

The cross section and its energy dependence is well described by STARLIGHT, while the GDL calcula¬ 
tion overpredicts the cross section by about a factor of 2. The agreement with STARLIGHT is somewhat 
surprising since its Glauber calculation does not include the elastic part of the total cross section, which 
is included in the GDL model. It has been argued that coherent p° production off heavy nuclei may 
probe the onset of the Black Body Limit, in which the total p°-nucleus cross section approaches 2kR\ 
at high energies [15]. The results from STAR and ALICE do not favour this picture. The cross section 
is instead reduced by about a factor of 2 compared with the GDF model [35], independent of energy, 
indicating that further work is needed to understand this process. It should be noted that none of the 
models in Fig. 5 include cross terms such as p + N —>■ p' + N. 

The photonuclear cross section, a(y+Pb —> p 0 + Pb), in STARFIGHT is almost energy independent for 
WyN >10 GeV. The increase in the Pb-Pb cross section, a(Pb + Pb —> Pb + Pb + p°), with ^Ann is thus 
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Selection Number of events Fraction STARLIGHT GDL 


All events 

7293 

100% 



OnOn 

6175 

84.7±0.4(stat.)i'j ) ^(syst.) % 

79 % 

80% 

Xn 

1174 

16.1±0.4(stat.)^Q5(syst.) % 

21 % 

20% 

OnXn 

958 

13.1±0.4(stat.)+°|(syst.) % 

16 % 

15 % 

XnXn 

231 

3.2±0.2(stat.)^gj(syst.) % 

5.2% 

4.5 % 


Table 2: The number of events that satisfy various selections on the number of neutrons detected in the ZDCs. 
OnOn corresponds to no neutrons emitted in any direction; Xn to at least one neutron emitted in any direction; OnXn 
to no neutrons in one direction and at least one neutron in the other direction; XnXn to at least one neutron in both 
directions. For the relative yield the systematic error is estimated, as explained in the text. 

almost entirely due to the increase in the photon flux at higher collision energies. 

The model by GM, although in agreement with the current result, has been criticized for using the Color 
Dipole model for a soft probe like the p° [35]. A recent publication shows that the calculation indeed 
has large uncertainties arising from the choice of p° wave function and dipole cross section [37]. 

The number of events satisfying the different fragmentation scenarios as well as the ratio to the to¬ 
tal number of events are shown in Table 2. The table also shows the expected fractions from the 
STARLIGHT [38] and GDL [35] models. These models assume that the probabilities for exchange 
of multiple photons in a single event factorize in impact parameter space. One should note that some 
of the fractions are correlated: the sum of (OnOn) and (Xn) should be 100%, and the sum of (OnXn) 
and (XnXn) should be equal to (Xn). This is the case within errors, but the sum is not exact, since the 
incoherent contribution is subtracted for each selection separately. The results in Table 2 are consistent 
with both the STARLIGHT and GDL calculations within three standard deviations. 

6 Conclusions 

The first LHC measurement on coherent photoproduction of p° in Pb-Pb collisions at ^Ann = 2.76 TeV 
has been presented. Comparisons with model calculations show that the measured cross section is in 
agreement with the predictions by STARLIGHT [16] and Goncalves and Machado (GM) [34], de¬ 
spite the idiosyncrasies in these models mentioned above. The Glauber-Donnachie-Landshoff (GDL) 
model [15, 35] overpredicts the cross section by about a factor of two. Comparisons with results from 
Au-Au collisions at RHIC energies indicate that this factor of two difference is independent of colli¬ 
sion energy in the range ^/vnn = 62.4-2760 GeV. In a recent preprint, it is argued that inelastic nuclear 
shadowing combined with the inclusion of intermediate states with higher mass in the y-vector meson 
transtition could explain the discrepancy [39]. Regardless of whether this is the correct explanation or 
not, it indicates that non-trivial corrections to the p° photoproduction cross section may become impor¬ 
tant at high photon energies. 

The relative yields for different fragmentation scenarios are found to be in agreement with predictions 
from the STARLIGHT and GDL models. This is important not only to confirm the assumptions in the 
two models but also because some experiments, e.g. PHENIX [40], have relied on a ZDC signal to 
trigger on ultra-peripheral collisions. To be able to relate such measurements to a photonuclear cross 
section, it is imperative that the probabilities for exchange of multiple photons are well understood. 

The total cross section is found to be about half the total hadronic inelastic cross section. This is an 
increase of about a factor of 5 from Au-Au collisions at Asnn = 200 GeV, where the fraction was about 
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10%. If the increase of the coherent p° photoproduction cross section continues to follow STARLIGHT, 
one can expect it to exceed the total hadronic production cross section of heavy ions such as lead or gold 
at a y^NN of about 20 TeV. 
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